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We propose a lattice gas model to account for linear chain structures adsorbed on (112) faces of 
tungsten and molybdenum. This model includes a dipole-dipole interaction as well as a long-ran- 
ge indirect (oscillatory) interaction of the form ~ cos{2kFr -\- •~p)lr, where kp is the wavevector of 
electrons at the corresponding Fermi surface and 95 is a phase shift. It is assumed that the structures 
are stabilized by an attractive indirect interaction along the chains. 

We have explicitly demonstrated that the periodic ground states strongly depend on a competition 
between the dipole-dipole and long-range indirect interactions. 

The effect of temperature in our model of linear chain structures is studied within the molecular- 
field approximation. The numerical results clearly show that for the dipole-dipole interaction only , 
all long-periodic linear chain phases are suppressed to low temperatures while phases with periods 
2, 3, and 4 dominate the phase diagram. However, when the long-range indirect interaction becomes 
important, the long-periodic linear chain phases start to fill up the phase diagram and develop a 
high thermal stability. 

We have chosen model parameters in order to reconstruct a sequence of long-periodic phases (for 
coverages less than 0.5) as observed experimentally at T = UK for Li/Mo(112) and Li/W(112). 
It would be interesting to verify our model and assumptions by checking experimentally the corre- 
sponding phase diagrams. 
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I. INTRODUCTION 

Chemisorption on metal surfaces has been attracting a 
great deal of attention for over three decades (for a review 
see, for example, [H and references cited therein). 
Its nature, however, despite quite involved experimental 
techniques and a number of data as well as many theo- 
retical attempts is far from being well understood. One 
of the recent challenging problem has been concerned 
with structures and phase transitions in metal submono- 
layers adsorbed on (112) faces of tungsten and molyb- 
denum ||l|,^. In particular, it was found from LEED 
experiments that many alkaline, alkaline-earth, and ra- 
re-earth elements adsorbed on these substrates form or- 
dered structures which, for low coverages, consist of lin- 
ear chains of adatoms being perpendicular to the atomic 
troughs at the surface. The chains of adatoms are of- 
ten far apart from one another in the direction along 
the atomic troughs (up to 9 lattice constants of the sub- 
strate) thus forming long-periodic (linear) chain struc- 
tures. Moreover, it turns out that the structures possess 
a high thermal stability Q . 

At higher coverages, however, submonolayer structures 
become much more complex for, in addition to coher- 



ent structures, there are one-dimensional incoherent or 
stressed submonolayer structures [Q. 

It is clear, that lateral interactions of adatoms are 
most important in determining the properties of ad- 
sorbed overlayers like their structures, phase transitions, 
thermal stability, etc. |jl[. These interactions at (112) 
surfaces of bcc metals (W, Mo) are believed to be due to 
a repulsive (isotropic) dipole-dipole interaction between 
adatoms with dipole moments, as well as to a long-ran- 
ge indirect (anisotropic) interaction via substrate elec- 
trons j^,^ . The delicate balance between the two interac- 
tions seems to account, at least for low coverages, for the 
observed properties of adsorbed linear chain submono- 
layers. So far, only qualitative (geometrical) arguments 
have been put forward in an attempt to understand to 
what extent the indirect electronic interaction might be 
responsible for the observed chain structures 

The purpose of this paper is to study how a compe- 
tition between the (repulsive) dipole-dipole and (oscilla- 
tory) long-range indirect electronic interactions could in- 
fluence ordered submonolayer structures as well as their 
equilibrium thermodynamic properties (phase diagram). 
This is done within a plausible two-dimensional lattice 
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gas model which might describe hnear chain structures 
for coverages up to 0.5 . 

The model is described in detail in section ||. In 



section [11 we calculate numerically the ground states 
(T = OK) and discuss their dependence on the model 
parameters. The mean-field approximation is used in sec- 
tion |^ to calculate the corresponding phase diagrams. A 
possible application to describe low coverage linear chain 
submonolayers of Li absorbed o n W( 112) and Mo(112) 
has been presented in subsection IV B . Finally, section ^ 
contains a general discussion and conclusions. 



II. THE MODEL 



It is generally believed, that chemisorbed atoms oc- 
cupy, at low coverage at least, only preferred sites of a 
lattice commensurate with the underlying substrate lat- 
tice. Here, we consider the (112) surface of bee metals 
(W or Mo) and the corresponding lattice of adsorption 
sites g (see fig. 0). 



(a) 



(b) 
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FIG. 1. (a) A schematic view of the (112) surface of bcc 
metals, (b) The corresponding (rectangular) lattice of ad- 
sorption sites shown by the crosses (a — 3.16A for W and 
a = 3.UA for Mo) 



It turns out, that lattice gas models are quite useful in 
studying overlayer structures and their properties (e.g. 
refs [^,H). As usual, a lattice gas model introduces the 
occupation variable np at each adsorption site r: np = 1 
if the rth site is occupied by an adatom and np — 0, if 
not. Then, within the grand canonical ensemble, the 
Hamiltonian is defined as 



n 



(1) 



where L denotes the rectangular lattice of NxN adsorp- 
tion sites with periodic boundary conditions and fx is the 
chemical potential (shifted, here, by the binding energy 
of an adatom to the substrate). Of course, ^ controls the 
adatom coverage 



(2) 



with (...) meaning a thermodynamic average. We are 
interested in a lattice gas model describing linear chain 
structures as observed experimentally for low coverages 
(9 < 0.5) and, herefore, we can neglect many-body inter- 
actions. Since, even in this case, the reliable first-princip- 
le calculation of effective lateral interactions is very dif- 
ficult, following a detailed discussion in ref. we have 
considered a pairwise adsorbate interaction F(r^ consist- 
ing of electrostatic and indirect interactions. Thus, we 
have assumed 



2(P 

V(r) = — 



Acos{2kF\y\+ip)^^S{x,0) for y ^ 



EbS{\x\,ai) 



for y = 
(3) 



where r = (x, y) and S stands for Kronecker's symbol 
(for notation see also fig. 1). These interactions are of es- 
sential importance in forming linear chain structures of 
alkalis, alkaline-earths, and rare-earths adsorbed on the 
(112) surface of W or Mo and they have the following 
meaning: 

(i) The first component in eq. (^, 2(P/\r\'^ , describes 
a repulsive dipole-dipole interaction between two 
identical adatoms adsorbed at r and r', respectively 
(see, for example, refs g and ||]). The dipole mo- 
ment of an adsorption bond, d, is usually estimated 
from the Helmholtz formula. Here, we neglect all 
depolarization effects due to many-body interac- 
tions and this seems to be justified for 9 < 0.5. 



(ii) The following component for y 7^ 0, Acos{2kp\y\ + 
ip)-^^5{x,0) , represents an asymptotic part of the 
indirect interaction between adatoms via conduc- 
tion electrons of the substrate (all relevant liter- 
ature on this interaction could be found in refs 
and ||l^). This particular interaction is highly 
anisotropic and is closely related to the existence of 
nearly flattened segments of the Fermi surface be- 
ing perpendicular to the [ 111 ] axis directed along 
the atomic troughs of the substrate. An amplitude 
A and a phase ip can be treated as phenomeno- 
logical parameters and kp denotes a wavevector of 
electrons at the Fermi surface. In our model we 
assume, following refs [|l]jl||, that kp = OAlA''^ 
{kp — 0.47A~^) for the tungsten (molybdenum) 
substrate. It is needless to say, that this indirect 
(oscillatory) interaction is long-ranged (^ 
and it could influence, as usually expected, proper- 
ties of the absorbed submonolayers. We have ne- 
glected all contributions to the indirect interaction 
related to hole segments of the Fermi surface [|[^ . 
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(iii) Finally, for y = the Ei,5{\x\, ai) term has been 
introduced to eq.(^ to facilitate the formation of 
linear chains of adatoms along the a;-direction. This 
is a part of an attractive indirect interaction be- 
tween adatoms along a chain and we assume that 
E}, = — O.OSeV at the nearest-neighbour distance ai 

§■ 

The model interactions, eq. (3), requires, however, some 
justification. 

First, the indirect, substrate mediated interaction (ii) re- 
flects only the leading "effective" dimensionality of the 
electronic states (e.g. refs [0, Q- jl^) mediating the 
interaction in the direction perpendicular to the chains. 
Here, we neglect a usual part of the indirect interaction 
decaying like ^ and associated with the " spherical" seg- 
ments of the Fermi surface [plp^]. Of course, r denotes a 
distance between two adatoms. 

Secondly, we would like to comment on what role sur- 
face electronic states may play in the indirect interaction. 
It is known, that they could contribute to this interac- 
tion provided the Fermi energy level intersects a band 
of the surface electronic states. Then, depending on the 
effective dimensionality of the surface states, the indi- 
rect interaction becomes proportional either to ^ (quasi- 
one-dimensional surface states jH- J15|) or to -i^ 
(two-dimensional surface states [|15[- |[16[). A period 
of the oscillations will be determined by a wavevector 
at the Fermi surface corresponding to the surface elec- 
tronic states. Recently, the field and photofield emission 
experiments on the clean W(112) surface have revealed 
the existence of surface states approximately 0.3eV be- 
low the Fermi energy jl^]- To our knowledge, how- 
ever, there is so far no experimental evidence that sur- 
face states mediate the indirect interaction in the case 
of the linear chain submonolayers on the (112) surface 
of W(Mo). Consequently, we have neglected a would be 
contribution to eq. (3) from the surface states. 
Thirdly, the model interaction (iii) could be thought of 
as a usual substrate mediated interaction in the form 
^ coa(2kpx + (j))S{y,0)/\xf (see, for example, refs 
p9{- p^). In order to describe the linear chain structures 
we might assume that (j) causes this indirect interaction 
to be attractive at the nearest-neighbour distance ai and 
to decay rapidly with distance. We cannot also exclude 
a possibihty that (ui) might be due to a virtual elastic 
distortion of the substrate surface [|3| . 

At this stage it is important to notice, that although 
the interaction described by eq. (||) is long ranged, we 
may restrict ourselves to an extended but finite range of 
interaction, R. A more detailed discussion will be pre- 
sented in the following section. 



III. PROPERTIES OF THE GROUND STATES 

A search for all the lowest energy configurations at 
T = OK for the model described by eqs (|l|) and (||) 



presents a formidable problem. Even in a simpler case of 
the dipole-dipole interaction only, as far as we know, it 
has not been possible up to now to find all the ground 
states. (For a review of the existing theoretical and nu- 
merical results concerning the ground states in two-di- 
mensional lattices see, for example, refs ||2^- and the 
literature cited therein.). However, let us note that the 
ordered phases found experimentally at T = 77K, as we 
have already discussed, are often composed of parallel 
linear chains, at least for low coverages {9 < 0.5). This 
suggests, that a search for the minimal grand canoni- 
cal ensemble energy configurations (the ground states) 
could be performed in two steps. First, we find the effec- 
tive Hamiltonian describing interactions between linear 
chains and only then we look for the ground states. 



A. The effective Hamiltonian 

Making use of the periodicity along the chains, we can 
write down the Hamiltonian as an effective Hamiltonian 
of a one-dimensional lattice gas type model, which de- 
scribes the interaction between linear chains (simply re- 
ferred to as chains, hereafter). This Hamiltonian has the 
form 



(4) 



where aj denotes chain variable with aj = 1 or depend- 
ing on whether there is a chain on a site j along the y 
direction or not. The chain-chain interaction energy at 
a distance a2\l\ reads 



+ 00 



(4a) 



where 



Vd{k,l) 



2<P 



cos(2fci?a2|/| -I- tp) 



The binding energy per chain can be written as 
-C(3)+i?6 



(4b) 



Here, C(3) = 1.202056. .. is the Riemann zeta-fimction 
and j, 1, and k are integers. 

It is well known, that the lattice gas effective Hamilto- 
nian, Heff, is analogous to the Hamiltonian of the Ising 
model in an applied field. For a class of long-range in- 
teractions (positive, convex, etc.) it has been possible to 
find in a one-dimensional case all the ground states for 
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any rational < - < 1 {p, q are integers with no com- 
mon multipliers) ^]. Moreover, the (0, /i) phase diagram 
represents the complete devil's staircase with a rather 
involved fractal behaviour Unfortunately, the 

effective Hamiltonian, Ti.ee, contains the indirect (oscil- 
latory) interaction and, therefore, we are not able to find 
the ground states rigorously. Also, the results of refs 
pO| , pi| do not seem to be applicable because the consid- 
ered interaction is truly long-ranged. 



B. A search for the ground states 

As usual, the long-range nature of the effective interac- 
tion between chains described by eq. (|4^) poses a prob- 
lem. It is easy to see that the interaction is dominated 



at large distances by the (oscillatory) indirect term and 
one has to sum up the infinite series (see eq. (^) for any 
would-be minimal energy configuration of u's. Instead, 
we approximate the problem (see also refs |^ and [ ^ ) 
by assuming that there is a sufficiently large range of the 
adatom - adatom interaction, say, R a2, beyond which a 
configuration under study is replaced by the average one 
having the same coverage 9. This average chain structure 
consists of equidistant "chains" with a spacing ao = a2/0. 
Moreover, adatoms are uniformly distributed along these 
"chains" with a density 1/ai. 

Finally, we restrict ourselves to periodic configurations 
of (t's, i.e. there is a period p (positive integer) such that 
aj-^-p ~ Uj for any integer j. These are called p-periodic 
configurations of cr's. 

Now, the energy per chain of a given p-periodic config- 
uration CToj • ■ • J fp_i can be written in the form 



1 R 



1 

E[ao, . . . , (Jp-l] = - X! X! ^(Oc^jC'modp 



U+i) 



j=0 1 = 1 



P 



p-i 

- £c) 



Er{9) 



(5) 



where 



y(0 



k{i} 

E 

k=-k(l) 



(6) 



with k{l) = Int(f^Vi?^ - /2). 

Here, the remainder of the energy per chain reads 



Er{e) = 



where 



A 



lnt(B.0) 

E 

n=l 



1 _ IJl. 



Re) 



+ 61^— [/i(a)cos^ + /2(a)sin^ + r(a)] , (7) 



02 



p-1 

;E^j- 



The first two components in eq. (^ depend on 
CTo: ■ ■ • I Cp-i and involve exact summations. The remain- 
der of the energy per chain, however, depends only on 9 
as it should have been expected. It is interesting to note 
that this model removes the particle- hole symmetry |^ . 

In the following, we propose a numerical procedure to 
find the (periodic) ground states. First, by using the bit 
representation of integers from the interval [2^'"^, 2^ — 1] 
we generate numerically all p-periodic configurations of 
cr's with p = 1, . . . ,Pmax- The number of configurations 
could then be reduced by making use of the translational 
and/or inversion symmetries. Secondly, we check explic- 
itly which configurations of cr's afford the minimal value 
to the corresponding i?[cro, . . . , crp_i] (p = 1, . . . ,p„iax)- 
Of course, the results depend on model parameters, such 
as d. A, ip,kp,R and Pmax- This will be discussed in the 
following subsection. 



a 



2kpa2 



C. The results 



/i(a) = — ln(2sin— ) for a = mod27r(a) > , 



Ol — TT 

/2(a) = — - — for a — mod2,r(Q:) > , 



Int(i?,e) 
■^-^ COS(nQ! + G9 — tt) 

n 

n=l 



To describe a competition between the dipole-dipole 
and indirect interactions we require three independent 
parameters: ^,<P, and kp. This results from eq. ^ 
where E[cro, ■ ■ ■ , crp-i] and fi + ec are rescaled by a fac- 
tor of A^^, thus changing effectively only the stability 
intervals of the ground states but not the ground states 
themselves. Moreover, our approximation concerning the 
energy per chain of a given p-periodic configuration of 
a's (eq. (|^)) and the numerical procedure introduce two 
more parameters: R and Pmax- In principle, calculations 
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should be performed for very large R (the larger the bet- 
ter) and for Pmax going to infinity to allow for any peri- 
odic structure. In practice, we have carried out calcula- 
tions for i? = 150 and Pmax = 23. This is a compromise 
between the computing time and more precise results. 
Our numerical tests show that relative changes of energy 
at Ra2 — 150a2 (in comparison with Ra2 — 1600a2) are 
less than one percent for a vast majority of p-periodic 
configurations of a's {p = 1, . . . ,Pmax = 23). Also, the 
results for larger values of Pmax change the results only 
quantatively. We believe, that to understand the role 
the indirect interaction plays in determining the ground 
states, no additional significant insight can be achieved 
by extending the numerical computations to higher val- 
ues of Pmax ■ 

We have applied our numerical procedure to determine 
the ground states for a number of different parameters: 

,2 

^,</?, and kp. Some of the results are presented in ta- 
bles ||-|| and fig. 2. The tables contain sequences of 
ground states. Here we use standard notation for the 
ground state, i.e. q/p, which means q chains in a unit 
cell of p sites along the y direction (fig. 1). At the same 
time ^ = p Bind in the present analysis we are not in- 
terested in what are the actual structures of the ground 
states. Also, it turned out to be convenient to express 
the strength of the indirect interaction. A, in miits of 
-0.137eVA (see refs §] and @). 

In fig. 2 we have shown explicitly examples of stability 
intervals of the corresponding ground states with respect 
to a reduced chemical potential ij,/ij,o (mo = ksT^T — 
lOOif). Note that only halves of the stability intervals 
for 1/2 and 2/4 have been depicted. As usual, the cov- 
erage as a function of reduced chemical potential forms 
the so-called " staircase" . 

Now, we shall discuss the presented results paying a 
particular emphasis to the role which the indirect inter- 
action plays in determining the ground states. 



i) Dependence on ^ 

This parameter is a measure of a competition between 
the dipole-dipole and indirect interactions. The results 
shown in table | and fig. 2 indicate that by decreasing ^ 
(or by increasing the strength of the indirect interaction) 
one causes the following changes to the ground states: 

,2 

• For very large values of ^ we practically ob- 
tain the complete devil's staircase which is shown 
in fig. 2(a) (see also refs |26|- [|§). 

,2 

• For large values of ^ many narrow long-periodic 
ground states (most of which correspond to higher 
coverages close to 0.5) disappear from the staircase. 
The remaining ground states usually change their 
stability intervals and might transform into the 
ground states (with a close coverage) (see fig. 2(b)). 



• For smaller values of ^ the narrow long-periodic 
ground states (corresponding to low coverages) con- 
tinue to disappear and we are left with few ground 
states some of which might be the long periodic 
ones (see fig. 2(c)). There can even be one ground 
state (for example, 7/20 for d — 0.3, A — 1, and 
other parameters as in fig. 2). This qualitative ar- 
gument is supported by our calculations for other 
values of (f and kp. For example, it also applies 
to ip — Q and ip — I.Btt with kp = 0.4lA~^ (see 
table |) but ^ must be larger than 2.25. 

At this point we would like to note that in the case of 
a sufficiently strong indirect interaction for O.Ott < p < 
1.57r and 0.40A~^ < kp < 0.50A~^ there is only a coexis- 
tence of the low-density (6* = 0) and high-density (9 = 1) 
disordered ground states. This means that an attraction 
effectively prevails between chains, thus leading to forma- 
tion of two-dimensional islands. This particular result, 
however, seems not to be accounted for by our model of 
linear chains (see section ID). 



ii) Dependence on Lp and kp 



It is easy to sec from tables g-y that the ground states 
depend on p and kp in a crucial way. The phase, 
(f, simply shifts uniformly the locations of minima of 
the corresponding indirect interaction thus allowing for 
the ground states with possible new periodicities or/and 
changed stability intervals (in /i). The wavevector, kp^ 
distorts uniformly the minima locations and a formation 
of the ground states is in this case even more complex. 
Indeed, the results presented in table || indicate that for 
(/? = each sequence of the ground states consists of 
the ground states having one (two) long-range period(s). 
Within our model this observation could be understood 
qualitatively by noting, that for a given kp^ the actual 
ground state q/p wins the energy competition between 
two terms: E[aQ, . . . , Op-i] — Er{9) involving a structure 
of a would-be ground state (i.e. (Tq, • ■ • 7 cTp-i) and Er{9), 
where 9 = - V^~n is such that a is close to or 27r 
(see eqs (||), (^, and (0)). This argument seems to be 
vaHd also for p e (0,^] and ip € [^,27r) despite the 
fact that the ground states start to gain on their ener- 
gies. And this gain for 1^9 S (|, ^) is even more dramatic 
for now there are several new ground states. Apart from 
long-periodic ground states there are also such ones as, 
for example, 1/2, 1/3 , and 1/4 (see tabled). 



IV. PHASE DIAGRAMS. MEAN-FIELD 
TREATMENT 



As a first step we adopt the mean-field approximation 
(a discussion concerning the validity of this approach is 
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postponed until the following section) to describe (equi- 
librium) phase diagrams for linear chains structures with 
the competing dipole-dipole and long-range indirect (os- 
cillatory) interactions. In the following, we restrict our- 
selves to the effective Hamiltonian Heff, eqs (Q), (4a), 
and (4b) with a finite range of the adatom - adatom in- 
teraction Ref f ■ The reason is that thermal fluctuactions 
" screen out" the long-range indirect interaction and it 
is plausible to assume for a qualitative analysis, at least, 
that Ref f does not depend on temperature and the dipo- 
le-dipole interaction can be neglected for distances larger 
than Reffa2. 

Following the standard procedure (e.g. ref. Q|) we can 
write the mean-field Hamiltonian in the form 



n. 



MF 



where 



and 



(8) 



(8a) 



(8b) 



Again, Uj denotes the chain variable with <Tj = or 1. 
Moreover, V{i ~ j) is given by eq. (^ with R replaced 
by Reff and fl = n + Sc with ec defined via eq. (4b). 

First, we solve the self-consistent eqs ( |8a| ) and (pb|). 
This becomes possible to do, if we assume, that the re- 
sulting structure is a p-periodic one, p = 1, . . . ,pmax- 
Usually the solution will not be unique (i.e. several or- 
dered structures are possible). Secondly, having found 
the solution(s) we can calculate the corresponding grand 
canonical ensemble potential(s) per chain 



1 



1=0 



1 



:Vd<J,. 



(9) 



Hence, the stable structure corresponds to the min- 
imal value of fl(T,fl), while the other structures are 
metastable or unstable ones. 

In the computations we assume Reff — 36 (i.e. 
Reff 0-2 ~ lOOA) and Pmax = 12 according to the "lever 
rule" between a computing time and more refined phase 
diagrams. It seems that Pmax = 12 is sufficient to ac- 
count for some linear chain structures determined by the 
LEED technique B. 



A. Dipole-dipole interaction 

We start with a model of linear chain structures which 
order via the dipole-dipole interaction only (with an ad- 
ditional attraction along the chains). This model might 



be relevant to adsorption on furrowed surfaces, where the 
indirect interaction can be neglected (small A or lack of 
the corresponding flattened segments of the Fermi surface 
of a substrate). 

The results of the calculations are presented in fig. 3. 
This is a temperature (T) versus reduced chemical poten- 
tial ( — ,/zo = fcsT, T = lOOA') phase diagram, where we 
have denoted equilibrium (stable) phases by their ground 
states notation q/p. Note, that there is one exception to 
the rule, i.e. the higher-temperature case, 4/10, which 
has no the ground state equivalent. The phase diagram 
is shown for /i > — 0.72/io and this is due to our restric- 
tion to Pmax = 12. The most interesting result is the 
following: the thermal fluctuations suppress long peri- 
odic phases and most of the phase diagram is dominated 
by the phases 1/2, 1/3, and 1/4. This occurs despite the 
long-range nature of the dipole-dipole interaction. (Note 
that a similar observation was reported in ref. [3^ ] for a 
two-dimensional model with short-range interactions). 



B. Application to Li/W(112) and Li/Mo(112) 

We claim that our model of linear chain structures 
might qualitatively account for the adsorption of lithium 
atoms on the (112) surfaces of W and Mo |], Q- [||. 
The adsorption bonds are polarized and at a low cover- 
age, d = 1.7D and d = IAD for the W and Mo sub- 
strates, respectively. Also, we expect the indirect inter- 
action to be important (see section ^) . Finally, no " hard 
core" effects (a relatively small atomic radius of a lithium 
atom) make the model more realistic. 

We have found two sets of the model parameters which 
might be relevant to Li/Mo(112) and Li/W(112) and 
the quantative features of the corresponding mean-field 
phase diagrams are summarized in figs 4 and 5. Both 
these figures have some common features. That is, the 
horizontal hatching indicates coexistence regions between 
ordered (equihbrium) phases denoted as in fig. 3. These 
are the so-called "mixed-phases" regions. A dashed curve 
represents a second order phase transition curve between 
the low density {9 — 0) lattice gas phase and 1/2 phase. 
This curve meets, at some angle, two tangent first order 
phase transition curves at a multicritical point P. 

Now, a direct comparison of the phase transition 
schemes for Li/Mo(112) [|| and for Li/W(112) || with 
a low temperature vs coverage sequence of transitions of 
fig. 4 and fig. 5, respectively, shows a qualitative similar- 
ity. 

In fig. 4, a low temperature scheme of phase transi- 
tion corresponds quite well to Li/Mo(112) because the 
first order transitions occur through forming islands of a 
phase in the "mixed-phase" region. It is also interesting 
to note a high thermal stability of 1/4 and 1/2 phases at 
the corresponding stoichiometrics. 

The phase diagram presented in fig. 5, however, needs 
some more explanation. The cross-hatching depletes a 
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region where the results do not seem to be reliable due 
to our restriction to Pmax = 12. We have also found that 
at low temperatures and close to the disordered phase 
(9 < 0.05) our numerical iteration procedure of solving 
eqs ( |8a| ) and (^) turns out to be non-convergent (oscilla- 
tory behaviour). Therefore, we consider only {9 > 0.05). 
Even so, we recover (not too close to the disorder phase) 
a low temperature sequence of 1/4, 1/3, and 1/2 phases. 
Of course, the present theory cannot describe the phase 
transition between 1/4 and 1/3 phases occuring via a 
kind of statistical mixture of phases (cf. [^,0 ) . Let us 
note that below ^ 80K there are long-periodic structures 
(i.e. 1/6 and 1/9). Thus, providing thermal fluctations 
do not destroy the long-range order, we expect that these 
phases could be observed experimentally below 77 K. We 
are aware of the approximations we have made so far 
and, therefore, no more direct comparisons are going to 
be made. A discussion of the validity of this approach to 
the other linear chain structures will be presented in the 
following section. 



V. DISCUSSION AND CONCLUSIONS 

The main results obtained in this paper fall into two 
categories, i.e. 

1. Ground states analysis of the effective energy, eqs 
(||), (^, and (^ concerning the role of the long-ran- 
ge indirect interaction. The results were obtained 
numerically and this was done without any approxi- 
mations. We have shown that the interaction mod- 
ifies in an essential way the ground states as ob- 
tained from the dipole-dipole interaction only. The 
general tendency is such that by increasing an am- 
plitude of the indirect interaction, A, one reduces 
the number of ground states. The full description, 
however, requires other parameters of the interac- 
tion, i.e. kp and if. A detailed discussion is pre- 
sented in section [II C . 



2. Temperature versus coverage (or T vs fj.) phase di- 
agrams based on the effective Hamiltonian Tie//, 
eqs (^, (pa|), and (4b), and computed numerically 
within the molecular-field approximation. 

It seems that (ii) needs more discussion. It is well 
known that the mean-field treatment does not account 
for (thermal) fiuctuations correctly especially in low di- 
mensional systems [p]p^. As a consequence, the long- 
range order might be destroyed and a topology of phase 
diagrams and/or the order of various transitions are 
not often correct. This has been discussed and demon- 
strated explicitly for short-range interactions by com- 
paring mean-field results with the corresponding Monte 
Carlo simulations, real-space or finite-size renormaliza- 
tion group (cf. pq , ^ ). 

In our opinion, it would be very interesting to study 
the problem of fiuctuations in the case of long-range 



(also oscillatory) interactions in "2.5Dimensional" sys- 
tems |p5| , ^ . In this context, our mean-field results can 
be treated as zeroth-order approximation. 

The role of the indirect interaction has a great impact 
on the phase diagrams despite the fact that it is quite 
weak (^ ~ O.ObeV). Numerical computations show in an 
explicit way that this interaction is responsible for caus- 
ing longer periodic phases to develop high thermal stabil- 
ity. Figs 3-5 illustrate this situation. The above conjec- 
ture, although based on the mean-field theory, is corrob- 
orated by the recent Monte Carlo simulation Note, 
however, that the indirect interaction used in ref. p9| is 
characterized by only three parameters and has no- oscil- 
latory behaviour. 

We think, it is a challange to extend (if practically 
possible) Monte Carlo simulation or other more "exact" 
methods to long-range (oscillatory) interactions, which 
seems to be important in adsorption at crystal surfaces. 
Also, it would be very interesting to determine experi- 
mentally the whole (T, 9) phase diagrams to check theo- 
retical predictions. 

The present work was concerned with adsorption of 
lithium atoms on furrowed surfaces of Mo and W . How- 
ever, other alkaline (alkaline-earth or rare-earth) atoms 
are larger in diameter and, therefore, our theory will have 
to be modified for "hard-core" effects. A kind of modifi- 
cation is also needed to account for more complex struc- 
tures. 

This paper has dealt with the indirect interaction de- 
caying like i which is closely related to both the face 
and form of the Fermi surface of the substrate (W and 
Mo). In general, however, the indirect interaction has an 
asymptotic form ~ cos{2kpr -\- Lp)/r^ , with /3 € [1, 5] 
We would like to end with the following comment. It is 
now well known that the precise power law behaviour of 
this interaction might depend upon the form of the Fermi 
surface, surface electronic states, virtual elastic distor- 
tion of the substrate surface, etc. Therefore, it would 
be very interesting to study a competition between the 
dipole-dipole and indirect (/? > 1) interactions in view of 
possible physical applications. This problem turns out to 
be quite difficult for it requires a novel numerical method 
of convergence of series and the results will be published 
elsewhere. 
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FIG. 2. Coverage as a function of reduced chemical potential, ^ (/io = fcsT, T = IQQK), for ip = 1.27r, kp = 0.4lA ^ and 
a = 3.16A. (a) The devil's staircase: d = 1.5 and A = 0. (b) The devil's-like staircase: d = 1.5 and j4 = 1. (c) The "normal" 
staircase: d = 0.5 and A = 1. Here d and A are expressed in units of Debye and —0.137eVA, respectively. For a detailed 
discussion see the text. 



FIG. 3. Mean-field phase diagram: T versus ^ (fio = kBT,T = WOK). Here d = 1.5D, = -0.05eV, a = 3.141. b and c 
denote the corresponding insets. For a notation see the text. 



FIG. 4. Mean-field phase diagram (T versus e ) ior d = ID, A = -0.274eVA, kp = 0.471-1, ip = 0.93n, Eb = -0.05ey 
and a = 3.141. For an explanation cf. text. 



FIG. 5. Mean-field phase diagram (T versus 61) for d = 1.5£), A = -0.137eFl, kp = 0.411"!, ^ ^ i_26'k, Eb = -0.05ey 
and a = 3.161. For an explanation cf. text. 
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TABLE I. Sequences of the ground states as a function of phase, tp. Hero, kp = 0.4ll"^ (a = 3.16A) and kp = 0.47A 



(a = 3.I4A) are considered because of a possible application to Li/W(112) and Li/Mo(112). ( ^ is expressed in units of 
-7.29^D'^eV-^A-^). 
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TABLE II. Sequences of the ground states as a function of wavevector, kp, for = and a = 3.16A. For a discussion see 
the text. ( ^ is expressed in units of -7.299Z)^el/"^A"^). 
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